Polarization and Strong Infra-Red Activity in Compressed Solid Hydrogen 
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Under a pressure of ~ 150 GPa solid molecular hydrogen undergoes a phase transition accom- 
panied by a dramatic rise in infra-red absorption in the vibron frequency range. We use the Berry's 
phase approach to calculate the electric polarization in several candidate structures finding large, 
anisotropic dynamic charges and strongly IR-active vibron modes. The polarization is shown to be 
greatly affected by the overlap between the molecules in the crystal, so that the commonly used 
Clausius-Mossotti description in terms of polarizable, non-overlapping molecular charge densities is 
inadequate already at low pressures and even more so for the compressed solid. 



PACS numbers: 78.30 .-j, 62.50.+p, 77.84.-s, 71.15.Mb 



The stretching mode of an isolated hydrogen molecule 
does not absorb infra-red (IR) radiation, since the electric 
dipole moment vanishes by symmetry and remains zero 
upon stretching. In the molecular solid the symmetry is 
lower, and for some structures first-order IR activity of 
vibrational modes is allowed Indeed, some IR absorp- 
tion occurs in the broken- symmetry phase (BSP or phase 
II) [§. More surprising is the dramatic increase in the in- 
tensity of IR absorption which occurs when crossing the 
boundary between the BSP phase and phase III (H-A) 
at around 150 GPa At 85 K the vibron oscillator 

strength goes from / ~ 0.0002 at 140 GPa to / ~ 0.0134 
at 167 GPa §. 

A considerable effort has been made to try to explain 
this remarkable behaviour. Hemley and co-workers 0] 
proposed that the increase in molecular overlap leads in 
phase III to the formation of charge-transfer states be- 
tween neighboring molecules, with the coupling between 
a vibron and charge-transfer excitations giving rise to 
charge oscillations between the molecules. Baranowski 
H invoked a strong ionization of the molecules in phase 
III. Mazin el al |J neglected the molecular overlap and 
treated the molecules as point objects which become po- 
larized in the quadrupolar field of the other molecules 
(EQQ model Q ) . Edwards and Ashcroft [jll| suggested 
that a dielectric instability associated with a charge- 
density wave causes an enhancement of the molecular 
dipole moments in phase III. Some theories || focus on 
purely static charge, others (?]] on purely dynamic (i.e., 
displacement-induced |l2]| ) and the others [p|JTlll take the 
two contributions to be comparable. There is also no 
agreement on the relative importance of intramolecular 
H|]|ll]] and intermolecular [7] charge transfer || . To set- 
tle these issues, a careful treatment of the bulk electric 
polarization, which is the central quantity at play in the 
physics of IR absorption, is needed. 

For that purpose we have used the "modern theory 
of polarization" Ju|, which shows how to compute the 
polarization of a periodic insulating system as a Berry's 
phase derived from the electronic wave functions. The 



information about the polarization is therefore in the 
phases of the wave functions, not in the periodic den- 
sity alone, whose dipole moment is in general ill-defined, 
depending on the choice of unit cell |H. The density- 
based picture is valid only in the Clausius-Mossotti 
limit of localized, polarizable units with non-overlapping 
charge densities fl5| and breaks down whenever signifi- 
cant charge derealization prevents an unambiguous as- 
signment of the electron density to particular atoms, 
molecules or ions in the solid. Such overlap effects are 
known to play an important role in the absorption of light 
by lattice vibrations in ionic systems JTs| , which typically 
have a much larger overlap than molecular crystals |l2) . 
But at 150 GPa solid hydrogen has been compressed by 
almost a factor of ten Jl(|, hence the motivation for ap- 
plying the modern theory of polarization to this problem. 

Using the Berry's phase approach we have calculated 
the Born effective charge tensors Z* 12 Fifl. In order to 



compare with experiments done on polycrystalline sam- 
ples, we used the expression for the oscillator strength 
f(j) of the contribution of a TO phonon mode j to the 
dielectric function averaged over all directions (lj] , 
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where £ a (j) = ^EkE^I/j W^ftj). k indexes 
the atoms in the unit cell of volume f2, ep (k;j) are the 
normalized eigenvectors associated with the mode, uj(j) 
is its angular frequency and m p is the proton mass. 

The crystal structure in the high-pressure phases has 
not yet been determined experimentally |l6| . For clar- 
ity we have investigated two of the simplest candidate 
structures for the compressed solid which have IR-active 
vibron modes Q] and are also energetically favorable 
]l8| , p9[ . The parameters a (lattice constant), c/a, bond 
length and 9 of Cmc2i and C2/m (see Fig. 1) were opti- 
mized by minimizing the enthalpy at fixed pressure using 
a variable cell shape method pQ] , and they are collected 
in Table I for several densities pltl. In some cases the 



1 



centers of the molecules were allowed to move away from 
the hep sites by symmetrical amounts in the yz plane, 
which does not lower the symmetry while allowing for 
the manifestation of the instability reported in Ref . [[ll| . 
The calculations were done in the local density approxi- 
mation (LDA) to the density functional theory 0] , with 
clamped nuclei. 




z=c/4 



z=0 



c/4 



FIG. 1. The Cmc2\ structure viewed along the c-axis (left) 
and in the yz plane (right). The centers of the molecules 
lie on hep sites, and the molecules in the two sublattices are 
tilted away from the c-axis by opposite angles 9 and —9. The 
C2/m structure is identical except that the two molecules are 
tilted in the same direction by an angle 9. 



Structure P r s c/a rbond 



fi 



in— phase /out 



C2/m 115 1.52 1.588 1.456 69.5 0.169 
C2/mt 115 1.52 1.583 1.460 70.6 0.589 



Although Z* x is non-zero by symmetry, we found it 
to be very small before the off-site relaxation. The pro- 
nounced anisotropy of Z* implies large dynamic charges, 
since displacements of rigid ions (static charges) give rise 
to isotropic diagonal Z* tensors |Q. The smallness of 
Z* x relative to the components in the yz submatrix sug- 
gests that most of the charge transfer occurs in the planes 
of the molecules, since a displacement Sx of an atom in 
the x direction changes the distances to other atoms in 
that plane to second order only. This analysis estab- 
lishes that intramolecular and/or in-plane intermolecular 
charge transfer are the dominant effects, but it does not 
distinguish between the two |||] . The small displacement 
away from the hep sites has a large effect on Z* , reveal- 
ing a strongly nonlinear dependence of the polarization 
on the positions of the molecules. 

Such large effective charges give rise to very strong 
IR absorption by the vibron mode in which the two 
molecules in the primitive cell vibrate out-of-phase (the 
in-phase mode is Raman active but IR-inactive in this 
structure [Q). Before the relaxation, the oscillator 
strength at 115 GPa is about 13 times larger than the 
experimentally measured value in phase III at 167 GPa, 
and it increases by a factor of 3.5 upon relaxation (see 
phascT able I). Due to closure of the LDA band gap, this struc- 
ture was not studied at significantly higher pressures. 
In Cmc2\ the symmetry of the structure dictates the 



Cmc2i 13 2.0 1.576 1.445 54.0 10~ fa 

Cmc2 1 115 1.52 1.576 1.445 54.0 0.00079 

Cmc2\ 115 1.52 1.574 1.451 56.6 0.014 

Ctoc2i 152 1.47 1.571 1.443 55.2 0.0026 

Ctoc2i 167 1.45 1.570 1.443 55.6 0.004 

Cmc2 1 180 1.43 1.569 1.443 55.9 0.0059 



Phase III 167 



undetected 0.0134 



9 x 10 B following form for the effective charge tensors: 
0.0114 

a \ 
b c i 
Ode) 
' 



0.049 
0.026 
0.034 
0.044 



Z*(l) 



JT(3) 



' Centers of the molecules are off-site (see text). The off- 
site displacement in a.u. of the upper molecule in Fig. 1 is 
-0.01 (C2/m), and 8y = -0.07 (Cmc2i). 




Z*(2) 



Z*(4) 




Sy = -0.1, Sz = 



TABLE I. Optimized parameters (except for r 3 = 2.0) and vi- 
bron oscillator strengths for several pressures. Pressures are 
in GPa and are obtained using the equation of state of Rcf. 
|l6| , angles are in degrees, and the other quantities in atomic 
units (a.u.). The experimental results for phase III are taken 
from Ref. Q. "undetected" means that the Raman-active 
vibron is not observed in the IR spectra. 



2.0 

1.52 
1.52+ 
1.47 
1.45 
1.43 



-0.0178 
-0.022 
-0.050 
-0.020 
-0.018 
-0.017 



■0.0033 0.028 0.044 0.0005 0.042 0.040 
-0.14 0.069 0.193 -0.0088 0.208 0.271 



-0.30 
-0.20 
-0.229 
-0.254 



-0.014 0.091 
0.069 0.216 
0.068 0.223 
0.065 0.228 



0.0087 0.243 0.348 
-0.015 0.260 0.347 
-0.017 0.283 0.380 
-0.020 0.307 0.413 



' Centers of the molecules are off-site (see text). 

TABLE II. Effective charges in atomic units for the opti- 
mized Cmc2\ at several densities. 



The existence of a center of inversion between the two 
molecules in the unit cell of C2/m leads to the relation 
Z*(l) = -Z*{2) = -Z*(3) = Z*(A) between the effec- 
tive charge tensors of the four atoms. For the optimized 
structure at 115 GPa we obtained, for atom 2 before and 
after the relaxation away from the hep sites, respectively: 

/ < 0.005 \ / 0.05 \ 

Z* ~ 0.53 0.18 , Z* ~ 0.96 0.28 ] . 

V 0.08 0.07 / V 0.15 0.16 / 



The calculated values for several densities are shown 
in Table II. The anisotropy is again very large, which 
implies large dynamic charges, but the overall values are 
somewhat smaller than for C2/m. Both vibrons are IR- 
active in this structure Q : the out-of-phase vibron has a 
large oscillator strength in reasonable agreement with the 
measured values in phase III (see Table I), whereas the 
in-phase vibron is about an order of magnitude weaker, 
which would still be enough to be detected experimen- 
tally, if Cmc2i was the correct structure for that phase. 
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The intensity of the absorption grows rapidly with pres- 
sure, in agreement with the behavior in phase III |§||0. 
At 167 GPa we obtain for the Szigeti effective charge q* 
H[l5| of the out-of-phase vibron q* ~ 0.047e, comparable 
to the measured value of q* ~ 0.032e Q, and its pressure 
derivative is dq*/dP ~ 2.3 x 1CT 4 e/GPa, close to the 
experimental value of ~ 2.1 x 10~ 4 e/GPa The effect 
of off-site relaxation of the molecules is again significant, 
with the oscillator strengths of the vibrons increasing 18- 
fold (in-phase) and 4-fold (out-of-phase) at 115 GPa. 




Position along bond 



FIG. 2. (a) electron density n(r) in a.u. for Cmc2\ along a 
line joining atoms 1(3) and 2(4) of Fig. 1. (b) asymmetry 
along that line, defined as a(£) = [n(£) — n(— £)] /2no, where 
no is the maximum electronic density in the hydrogen atom 
@. For r s = 3.4 a(£) was multiplied by (3.4/2. 0) 4 for com- 
parison with r s = 2.0. The molecular centers are on hep sites, 
and therefore the intramolecular a(£) is very small Jn|, and 
yet the vibron oscillator strength for r s = 1.52 is quite large. 

More insight into the physics of polarization can be 
obtained by looking at the electron density. In Fig. 2a 
we show the density along the axis of a molecule in the 
Cmc2i structure for several pressures. The asymmetry 
induced by the neighboring molecules can be clearly seen 
in the higher pressure curve, and is more visible in Fig. 
26, where we have plotted the asymmetry function a(£) 
[jllj . If the molecules can be considered to be separate po- 
larizable objects, as in the Clausius-Mossotti model, then 
the asymmetry of the molecular density will determine 
the polarization. For example, in the widely used EQQ 
model P, ^0|p4|j25] | the molecular dipoles are induced by 
the electric field at the centers of the molecules due to 
the quadrupolc moments Q of the other molecules. In 
this model the intramolecular asymmetry scales as 1/r 4 
(assuming that Q and the molecular polarizability do not 
change much with pressure |2(|). Indeed this is the case 
at low and moderate pressures ( r s > 2.0), as can be 
seen in Fig. 2b from the agreement in the intramolecu- 
lar region between the curve for r s = 2.0 and the curve 
for r s = 3.4 scaled by (3.4/2.0) 4 . However, the scal- 
ing clearly does not hold up to the megabar range, since 



the intramolecular a(£) would have to be around 3 times 
larger for r s = 1.52 than for r s = 2.0, whereas it is found 
to be slightly smaller and decreasing with pressure in that 
range. Therefore in the megabar range the EQQ inter- 
actions do not account for the local field acting on the 
molecules, contrary to the model used in Ref. |l. 

In Ref. jll] a(£) in the intramolecular region was in- 
terpreted as an estimate of the dipole moment of the 
molecule in the solid, which leads to the prediction that 
the vibron oscillator strength should scale roughly as the 
square of its magnitude. Instead, we have found that al- 
though the intramolecular asymmetry is comparable for 
r s = 1.52 and r s = 2.0, the oscillator strengths of the 
vibrons are ~ 800 (in-phase) and ~ 1300 (out-of-phase) 
times larger for r s = 1.52 (see Table I). The origin of this 
surprising result can already be seen in Fig. 2, where it is 
clear that for both these densities <x(£) is large in the in- 
termolecular region, where the density is small, which is 
related to the fact that the outer regions of the molecules 
are most easily polarizable p7j . Thus, when considering 
polarization effects, molecular overlaps may be impor- 
tant, even if they are small, causing the dipole moments 
of the molecules to be ill-defined and compromising the 
Clausius-Mossotti picture. 

At low enough densities the Claussius-Mossotti model 
must be correct; however, to our knowledge there has 
never been a thorough investigation of how low the den- 
sity must be for this limit to apply. This is now possible 
using the modern theory of polarization p3]. Here we 
present results for the Cmc2\ structure as a function of 
volume. (We emphasize that this is not meant to repre- 
sent real hydrogen at low pressure, where the molecules 
are not oriented but instead behave as quantum rotors 
p4| , p5[ , but is sufficient to establish the desired points re- 
garding the sensitivity of the polarization to the molecu- 
lar overlap.) In this structure there is a net spontaneous 
polarization Po along the c-axis, and in Fig. 3 we plot 
Po as a function of r s , together with the value given by 
the EQQ model and the dipole moment for a particular 
choice of unit cell. The cell we used goes from z = to 
z = c in Fig. 1, which is the obvious choice in the low 
pressure limit since the electronic density will vanish at 
the boundaries and the dipoles will become well-defined 
in that limit. For r s > 1.52 atoms 1 and 3 in Fig. 1 have 
an excess of electron charge with respect to atoms 2 and 
4, as inferred from a(£). In a Clausius-Mossotti frame- 
work this would lead to a net polarization (sum of the 
molecular dipole moments) along minus z, and that is in- 
deed the result of the EQQ model (and also of the dipole 
moment for r s > 2.5), but between r s ~ 1.7 (~ 50 GPa) 
and r s ~ 3.0 (around atmospheric pressure) the Berry's 
phase polarization has the opposite sign. Most surprising 
is the fact that the corrections to the models are so large 
even at low pressures (large r s ). Since the three quan- 
tities in Fig. 3 start to converge to each other only for 
r s > 3.0, we conclude that the Clausius-Mossotti limit 
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is reached only at negative pressures. A possible expla- 
nation for this observation is that at zero pressure there 
must still be some overlap of the wavefunctions in order 
to balance the attractive interactions, and only at nega- 
tive pressure does this overlap become truly negligible. 




r s (a.u.) 

FIG. 3. Spontaneous polarization per unit cell as a function 
of r 3 for the /icp-centered Cmc2\. The structural parameters 
are the optimized parameters for r s — 1.52 (see Table I). The 
Berry's phase result is compared with the value given by the 
EQQ model, using the quadrupole moment and polarizability 
of the H2 molecule quoted in Ref. |^5j , and with the "dipole 
moment" (see text) of the charge density in a unit cell. 

In summary, we have presented the first ab-initio calcu- 
lations of the intensity of IR absorption in solid hydrogen, 
showing that the intense vibron IR activity found in some 
structures is caused mainly by large dynamic charges. 
The strong dependence of the effective charges on the 
symmetry and structural parameters offers a plausible 
explanation for the difference in IR absorption between 
the two high-pressure phases, and the study of other 
structures with the present approach may help identify 
the crystal structure of those phases. The physics of po- 
larization of solid hydrogen appears to be more subtle 
and less amenable to simple models than has been pre- 
viously assumed. This seems to be caused by the very 
delocalized nature of the induced dipoles, which extend 
out to the exponentially decaying tails of the molecular 
charge density | p7| , rendering the bulk polarization very 
sensitive to even small overlaps between the molecules. 
These findings may be relevant for other molecular solids, 
since the Clausius-Mossotti picture is assumed in most 
treatments of the IR absorption in such systems |28| . 
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